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ABSTRACT 

We present 5.5-7.6 /im spectra of comets 71P/Clark (2006 May 27.56 UT, 
Th = 1.57 AU pre-perihelion) and C/2004 Bl (LINEAR) (2005 October 15.22 UT, 
Th = 2.21 AU pre-perihelion and 2006 May 16.22 UT, = 2.06 AU post- 
perihelion) obtained with the Spitzer Space Telescope. The z/2 vibrational band 
of water is detected with a signal-to-noise ratio of 11 to 50. Fitting the spectra 
using a fluorescence model of water emission yields a water rotational temper- 
ature of < 18 K for 71P/Clark and ~ 14 ± 2 K (pre-perihelion) and 23 ± 4 K 
(post-perihelion) for C/2004 Bl (LINEAR). The water ortho-to-para ratio in 
C/2004 Bl (LINEAR) is measured to be 2.31 ± 0.18, which corresponds to a 
spin temperature of 26^2 K- Water production rates are derived. The agreement 
between the water model and the measurements is good, as previously found 
for Spitzer spectra of C/2003 K4 (LINEAR). The Spitzer spectra of these three 
comets do not show any evidence for emission from PAHs and carbonate min- 
erals, in contrast to results reported for comets 9P/Tempel 1 and C/1995 01 
(Hale-Bopp). 

Subject headings: Comets: individual (71P/Clark, C/2004 Bl LINEAR): in- 
frared: solar system 
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INTRODUCTION 



Production rates of volatiles are used to quantify cometary activity and to measure 
volatile abundances. Water is the dominant ice in cometary nuclei and its sublimation 
governs the activity of comets at heliocentric distances, Vh, <3 AU from the Sun. Water 
production rates, Q(H20), the rotational temperature, Trot(K), and the nuclear spin tem- 
perature inferred from the water ortho-to-para ratio (OPR) are physical parameters useful 
in a variety of studies related to cometary atmospheres and cometary physics. 

Unlike other constituents of cometary atmospheres, water vapor is difficult to observe di- 
rectly from the ground because of terrestrial atmospheric absorption. Non-resonance fluores- 
cence vibrational bands (hot bands), not absorbed by telluric H2O, can b e observed from the 



groun d in relatively active comets {Q{B.20) ^ a few 10^*^ molecules s~^; iDisanti &: Mumma 
20081 ). The detection of the H2O fundamental bands of vibration {ui, z/3 at 2.7 fim, 1/2 at 
6.3 /xm), and of H2O rotational lines in the submillimeter and far-infrared domains, requires 
space-based instrumentation. Detection of these signatures has been achieved in the coma 
of several comets using spectrographs on the VEGA, Deep Impact, and Rosetta spacecrafts, 
and instruments on the Kuiper Airborne Observatory, the Infrared Spac e Observatory (ISO) , 



the Submillimeter Wave Astronomy Satellite, and the Odin satellite (lA'Hearn et al.l 



Biver et al.l 120071 : iBockelee-Morvan et al.ll2004l and references therein; iGulkis et al.l 120071 ). 



2005 



Gehrz et al. 


2007; 


Werner et al. 


2004) 



observations of the z/2 vibrational band of water vapor at 6.3 yum with a spectral resolution 
of ~100. This band was detected with Spitzer in comet C/2003 K4 (LINEAR) with a high 
signal-to-noise ratio (^50), which enabled a detailed study of the band shape using a fluo- 
rescence excitation model, assessment of the radial gradient of emission with cometocentric 
radius, and the measu rement of the rotationa l temperature and ortho-to-para ratio (OPR) 



of the detected water (IWoodward et al.ll2007l ). 



An important by-pro duct of the analysis of the Spitzer spectra of comet C/2003 K4 
(LINEAR) performed by IWoodward et al.l (120071 ) was a search for emissions in excess to 
the water band in the 5.8-7.6 /xr n domain, which i s rich in spectral signatures of gaseous 
compounds and minerals. Indeed, iLisse et al.l (l2006l . 120071 ) reported emission from carbonate 
minerals at 6.5-7.2 /xm and organic (PAH) emission at 6.2 /xm in the Spitzer spectrum of 
comet 9P/Tempel after the collision with the De ep Impact impactor, a nd in comet C/1995 
01 (Hale-Bopp ) based on their re-anal ysis of the ICrovisier et al.l (119971 ) ISO spectrum. The 
methodology of iLisse et al.l (l2006l . 120071 ) does not rely on a search of individual signatures but 
from the decrease of the residual after a fit of the whole (5.2 to 40 /im for Spitzer spectra) 
observed spectrum by a modeled spectrum which can incorporate many constituents. A re- 
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analysis of the ISO data of comet Hale-Bop p bvlCrovisier &: Bockelee-MorvanI (120081 ) did not 
confirm the detection of PAHs reported by lLisse et aLl (120061 ). and suggests that carbonate 
signatures are at most marginal. After detailed modeling of the water z/2 emission, the spectra 
of C/2003 K4 (LINEAR) do not show any evidence for emission from PAHs and carbonate 
minerals (jWoodward et al.ll2007l ). Carbonates have been identified in very tiny quantities in 



the d ust grains collected in the coma of 8 IP/ Wild 2 by the 5'tord?/5t spacecraft (Wirick et al. 



20071 ) . PAHs have been also detected in Stardust samples (ISandford et al.l 120061). but a sub- 
stantial portion of thi s material may not be of cometary origin (jSandford fc Brownlee 2007 : 
Spencer &: Zard 120071 ). PAHs were tentatively identified in near-IR and near-UV cometary 



spect ra, but these detections are controversial (see the discussion in iBockelee-Morvan et al. 
2004h . 



Here, we present longslit Spitzer spectroscopic observations of comets 71P/Clark and 
C/2004 Bl (LINEAR). Observations in the 5.3-7.6 /im spectral region are presented in 
§ [21 The d ata are analyzed with a H2O fluorescence model using the same methodology 
described in lWoodward et al.l (120071 ) and the results are presented in §[3l Residual emissions 
are discussed in § 13.61 We conclude in § HI 



2. OBSERVATIONS AND REDUCTION 



Comet 71P/Clark was discovered on 1973 June 09 by Michael Clark at Mount John 
University Observatory in New Zealand. This comet, with an orbital period of 5.52 years, 
belongs to the dynamical class of Jupiter-family comets. Its last perihelion approach was on 
2006 June 06.80 (perihelion distance q = 1.562 AU). Little information is available on this 
comet in the literature. The radius of its nucleus has been estimated to 0.68 ± 0.07 km from 
observations with the Hubble Space Telescope obtained at Vh = 2.62 AU, to 1.3 km from 
observations at r/j = 4.67 AU performed with the Keck telescope, suggesting an elongated 
nucleus with an axis ratio a/b ^ 2.85 (ILamy et al.l l2004l : iMeech et al.l |2004| ) . The orbit 
of this comet is peculiar. 71P has a non-zero, non-gravitational thrust in the direction 
normal to the orbital plan e (the A3 non-gravitational pa rameter) that is well-determined 
with good confidence (e.g., Sekanina 1993 : Nakano 2006a| jbl). Erratic time variability of the 
non-gravitational parameters over nucleus revolutions were interpreted as due to a change 
in the surface distribu tion of active areas and of the orientation of the spin axis (jSekanina 
I993I : ISzutowicg|l999l ). The water production rate of 71P/Clark at perihelion in 1973 and 
1984 was estimated to be 8-9 xlO^'' molecules s^^ based on its visual brightness (jSekanina 
I993I ). IWoodward et al.l (120081 ) find strong evidence for CO and CO2 gaseous emission near 
4.5 /xm in the extended coma pre-perihelion (r/j = 1.607 AU) during the 2006 apparition. The 
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Spitzer IRS observations of 71P/Clark were conducted ~ 10 days before its 2006 perihelion 
passage. 

Comet C/2004 Bl (LINEAR), first identified as an asteroidal object w ith a vi s ual ra ag- 
nitude of 19.1 on 2004 Jan 29.16 UT at a r/, ~ 8.1 AU, was discovered bv lYoungl kooi \ to 
be a comet on 2004 Jan 30.1 UT with a 3" coma, and an extended tail 4", PA ~ 310°). 



Orbital elements derived from astrometric observations ( IMarsdenl I2004J ) and subsequently 
refined over a data arc spanning years 2004 through 2007 (1779 observations), indicate that 
this comet is dynamically new, highly inclined to the ecliptic {i ~ 114°), with a perihe- 
lion of 2006 Feb 07.89 UT with a periapsis, q = 1.602 AU. Although a rela tively bright 



comet, few ground-based observations are currently reported in the literature. iBonev et al. 



(120071 ) conducted i?— band imaging observations of C/2004 Bl (LINEAR) over a ~ 60 day 
period between 2006 June and 2006 August, reporting that the come t underwent episode s 
of variable activity on timescales of days, with a dust production rate ( A'Hearn et al.lll984j ). 
Afp ~ 200 — 350 cm. The comet also exhibited a substantial coma 2 x 10^ km diam- 
eter) near the end of 2006 June, and slowly faded in i?— brightness, from ~ 14.3 mag to 
~ 15.1 mag, 193 days after perihelion. 



2.1. Spitzer IRS 

Comet 71P/Clark was observed once = 1.565 AU, a Spitzer-comet distance of 
0.920 AU, and a phase angle of 37.8°) in the second order of the short-wavelength, low- 
resolution module (SL2) on 2006 May 27.56 UT as part of a large Spitzer Cycle 2 comet 
survey, program identification (PID) 20021 (PI: C.E. Woodward), astronomical observation 
request (AOR) key 13818368, and processed with IRS reduction pipeline S15.3.0. The SL2 
slit is 3.7" wide and provides 57" of spatially resolved spectra (1.85" pixel"^) with a spectral 
dispersion of 0.06 /im. Nine spectra (60 s x 5 cycles) at 5.2-7.6 fim were recorded in a 3 x 3 
spectral map (executed with no peak-up), with 2.47" x 38.0" steps (perpendicular x parallel 
to the long slit dimension). The comet was also observed with the long- wavelength, low- 
resolution (LL) module of IRS, with which we used to verify the pointing of the spacecraft. 
Further analysis of the LL spectrum is beyond the scope of this paper. 

Comet C/2004 Bl (LINEAR) was observed twice, both pre- and post-perihelion with the 
SL2 module as part of a multi-cycle Spitzer initiative to investigate the heliocentric evolution 
and activity of a given comet, PID 20104 (PI: D.H. Wooden). The SL2 spectra discussed here 
were derived from IRS spectral mapping observations obtained on 2005 October 15.22 UT 
(pre-perihelion, Vh = 2.210 AU) AOR key 15787008, and on 2006 May 16.24 UT (post- 
perihelion, = 2.058 AU), AOR key 15791616. Pre-perihelion, three spectra (240 s x 5 
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cycles) were obtained in from a 3 x 1 spectral map, while post-perihelion five spectra (60 s 
X 3 cycles) were generated from a 5 x 1 spectral map. Both AORs were executed using a 
mapping footprint of 2.50" steps (perpendicular to the long slit dimension) with no peak- 
up on the comet nucleus, accompanied by appropriate shadow (background) observations. 
Spectral observations of comet C/2004 Bl (LINEAR) using the short-high (SH, 9.9-19.6 /im) 
module were also obtained, although detailed analysis of these spectra are beyond the scope 
of this manuscript. 

A reconstruction of the pointing for the 71P/Clark SL2 spectral maps indicates that 
there is an 8" perpendicular offset between the SL2 extraction position and the JPL Horizon^ 
ephemeris position of the comet nucleus at the epoch of our observations. Fig. [l](c). This 
offset is confirmed by comparing the peak in the observed surface brightness distribution in 
the long-low (LL2, 14-20 /im) cube which is nearly coincident with the Horizons position. 
Fig. [U^d). The spectral maps for both IRS pointings of C/2004 Bl (LINEAR) encompassed 
the peak in the surface brightness of the coma (Fig. [1]), hence SL2 extraction apertures 
centered on the nucleus are possible. 

The reduction of the spectr a and wavelength ca libra tion of the SL2 d a ta fol lows the 
formalism discussed in detail by iKelley et al.l (120061 ) and IWoodward et al. J2007h . "Hot- 



pixels" occasionally present in individual extractions derived from the spectral maps were 
removed and replaced by nearest neighbor interpolation, or ignored altogether. A summary 
of comet physical parameters, spectral extraction apertures, and derived water production 
rates (see § [3]) are given in Table [TJ 



3. DATA ANALYSIS 

3.1. Model Fitting of the U2 Water Band 

As for C/2003 K4 (LINEAR), features observed between 5.8 and 7 fim above the con- 
tinuum background in the 71P/Clark and C/2004 Bl (LINEAR) Spitzer spectra (Fig. W) 
coincide with the ro-vibrational structure of the V2 water band at the resolution of the short- 
low (SL2) spectrometer. In order to study whether water band emission fully accounts for 



1 in 


Woodward et al. 


( 


2007) 


Woodward et al. 


(2007 


)• 



^http: / /ssd.jpl. nasa.gov/horizons.cgi 
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In summary, we use the model of fluorescence water emission presented by lBockelee-Morvan fc Crovisie 



( 119891 ) which considers five excited vibrational states and their subsequent radiative cascades. 
Indeed, the V2 band, which dominates the water spectrum in the 5.8-7.1 /xm spectral region, 
is significantly populated by radiative decay from higher excited vibrational states. The 
resulting emission rate of z/2 is 2.41 x 10~^ s~^ at r/^ = 1 AU from the Sun. For computing 
the line-by-line fluorescence, 32 ortho and 32 para rotational levels are considered in each 
vibrational state. The rotational populations in the ground vibrational state is described 
by a Boltzmann distribution at a temperature Tj-ot- The gas expansion velocity, v^xp^ was 
fixed to 0.8 km s^^ and the water photodissociation rate was taken equal to 1.6 x 10^^ s^^ 
{r^ = 1 AU). These parameters influence to a small extent the retrieved water production 
rates. 

From the model output, synthetic Spitzer spectra are generated by convolving the in- 
tensity of the individual ro-vibrational lines with the instrumental spectral response of the 
spectrometer. The spectral resolution is determined to be 0.0614 ± 0.0027 /im from the 
pre-perihelion spectrum of C/2004 Bl (LINEAR), in good agreement with the SL2 spectral 
resolution (0.0605 /xm) cited in the Spitzer IRS Data Hand Book v. 3. 20. Our derived value 
is used for the analysis of the comet spectra in this paper. A slightly lower value of the SL2 
spectral resolution, 0.065 /im (determined from the wavelength calibration at the SL2 slit 
edge), was used in the analysis of the C/2003 K4 (LINEAR) spectra as sp ectral extractions 



of th is comet were performed near the edge of the slit (see § 2.2 and 3.1 in lWoodward et al. 



2003). 



The observed spectra are fitted with a composite curve consisting of the modeled water 
spectrum superimposed on a 5*^-degree polynomial. This latter polynomial determines the 
underlying continuum due to dust emission. The best-fit modeled spectrum is searched for by 
applying a least-squares method that uses the gradient- search algorithm of Marquardt. The 
only free parameters of the model are the water production rate Q(H20) (or equivalently, 
the water column density). Trot, the ortho-to-para ratio (OPR), and the polynomial for 
the continuum. However, the signal-to-noise ratio in the spectra of comets 71P/Clark and 
C/2004 Bl (LINEAR) post-perihelion is not high enough to constrain the OPR, so this 
parameter was fixed when analyzing these data. 

Figure [2] shows the model fits superimposed on the observed spectra. In the right 
panels, B, D and F, the continuum has been subtracted. Model fits with residuals are shown 
in Figs. Sand U for 71P/Clark and C/2004 Bl (LINEAR), respectively. The intensity of 
the water band in continuum-subtracted spectra and the retrieved Trot, OPR and Q(H20) 



^ |http : //ssc ■ spitzer . caltech.edu/IRS/dh 
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are given in Table [H 



As discussed in IWoodward et al.l (120071 ). the spectral resolution of Spitzer SL2 is to low 
to resolve the individual ortho and para water lines. Thus a fully independent estimate of 
the OPR and Trot cannot be derived. However, a robust determination of these parame- 
ters is possible from model fitting and minimization, because several para lines (near 
6.12 and 6.4 /xm) are well separated in wavelength from strong ortho lines, while several 
peaks, dominated by ortho lines, have their relative intensities essentially influenced by Trot- 
Spectra obtained with high signal-to-noise ratios are of course requisite. Figure M^a) shows 
iso-x^ contours in the (Trot, OPR) region that provides the best fit to the C/2004 Bl (LIN- 
EAR) pre-perihelion spectrum. The coefficients of the polynomial describing the continuum 
background have been set to the values of the best fit. The reduced obtained for the 
best fit is 0.97 in the 5.85-7 /xm region. The weak inclination of the ellipses (i.e., their 
roundness at the scales of the plot) indicates that the two parameters are poorly correlated 



(IBevington fc Robinsonlll992l ). The linear correlation coefficient is equal to 0.034. In Fig- 
ure[3]^b), red iso-x^ contours are computed by fixing Trot and OPR and then minimizing 
with the coefficients of the polynomial and the water production rate left as free parameters. 
With respect to the iso-x^ contours computed with the previous method (superimposed in 
black), these iso-contours extend a bit farther in Trot, while they coincide along the OPR 
dimension. This demonstrates that Trot is weakly correlated with the coefficients of the poly- 
nomial, whereas OPR is not correlated. Quoted uncertainties in Table [1] and throughout the 
paper are l-cx (68% confidence level) uncertainties computed from the diagonal elements of 
the covariance matrix. They correspond approximatively to the Ax^ = 1 red contour in 
Figur e [Mb). The joint confidence region for (Trot, OPR), delimited by the Ax^ = 2.3 con- 



tour (IBevington fc RobinsonI 119921 ). is not used here as these parameters are not physically 
linked. 



3.2. Water Modeling Results 

In the case of 71P/Clark, model fitting was performed for two different extraction aper- 
tures along the slit, 3.7"x 5.6" (top of Fig. iD and 3.7"x 17.9" (bottom of Fig. SD. For 
the small extraction aperture, there is significant disagreement at 6.35-6.40 fim (this region 
was masked during the fitting procedure to obtain successful results in the other domains). 
The intensity of the water band is less intense than expected from modeling, resulting in 
negative residuals. This low intensity is puzzling. At low temperatures, the 6.40 /im peak 
arises mainly from the Iqi — lio and In — 2o2 lines, whereas the 6.50 /im peak comes mainly 
from the Iqi — 2i2 and 2i2 — 22i lines. The ratio of their intensities /(6.40 /im)//(6.50 /xm) 
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decreases with decreasing Trot (see Fig. 5 of IWoodward et al.l 120071 ). but reaches a constant 
value of ~ 0.4-0.5 at low T and fluorescence equilibrium. For the large extraction aperture, 
the discrepancy is shifted towards the 6.45-6.50 /im domain (Fig. H]). Therefore it is likely 
that the misfit in the 6.35-6.50 /xm range is related to incorrect background subtraction 
and/or instrumental artifacts. Outside this spectral domain, the agreement between our 
model of H2O emission and the Spitzer spectra of 71P/Clark is good. 

The best fits to the spectra of 71P/Clark correspond to rotational temperatures Trot of 
5±13 K (large extraction aperture) and 14 ±4 K (small extraction aperture), both consistent 
with Trot < 18 K (similar values are found for OPR = 2.5 and 3.0). Indeed, the 6.50 /xm peak 



is les s intense than the 6.65 /xm peak only for Trot < 20 K (see Fig. 5 of IWoodward et al. 



20071 . where the emergent intensity of the H2O z/2 band is plotted as a function of rotational 
temperature). The reverse was observed in Spitzer spectra of C/ 2003 K4 (LINEAR) from 
which Trot values ~ 30 K were determined (IWoodward et al.l 120071 ). 



The signal-to- noise ratio achieved in the C/2004 Bl (LINEAR) pre-perihelion water 
spectrum (~50 in-band intensity on the 3.7" x 17.9" extraction aperture considered for model 
fitting) permits a determination of the water ortho-to-para ratio (OPR), as done from Spitzer 
spectra of C/2003 K4 (LINEAR) JWoodward et al.ll2007h . The best fit, shown in Figure [3 
corresponds to Trot — 14 ± 2 K and OPR = 2.31 ± 0.18. Also shown in Fig. [5] is the best fit 
obtained with the OPR fixed at a value of 3. The is increased by Ax^ = 10 with respect 
to the best fit obtained with OPR = 2.31. 

The agreement between our model and the Spitzer spectrum is excellent, except for 
the structure at 6.85 fim where the predicted intensity is lower than observed with a 2-a 
discrepancy. Taking the 6.05 /im peak as a re ference since its intens ity is weakly depen- 
dent upon Trot and OPR (cf.. Figs. 5 and 6 of IWoodward et al.l 120071 ) . the intensity ratio 
/(6.85 /im) //(6.05 /xm) slowly inc reases with increasing Trot (and decreasing OPR) (Figs. 
5 and 6 of IWoodward et al.l 120071 ). The observed ratio would be typical of Trot ~ 50 K. 
However, the relative intensities of the other peaks (in particular, 7(6.65 yum)/J(6.05 fim) 
and 1(6.18 /im)//(6.05 /xm)) are indicative of a low Trot- Model fitting to the 5.9-6.7 /xm 
part of the s pectrum yields the sarn e results as those obtained from the whole spectrum. As 
discussed in lWoodward et al.l (120071 ). the band regions most sensitive to the OPR lie at 6.12 
and 6.4 /im where the emission is dominated by para lines: model fitting with OPR = 3 
results in excess emissions centered at 6.12 /zm and 6.4 /im in the residuals. The discrepancy 
at 6.85 /im is likely related to instrumental effects or flaws in the background subtraction 
(the negative feature at 7 /im in the residuals is not explained either). 

The best-fit modeled spectrum for C/2004 Bl (LINEAR) post-perihelion is shown in 
Figure [51 For this model fitting, we fixed the ortho-to-para ratio to the value OPR = 2.3 
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derived from the pre-perihelion data. The derived rotational temperature is Trot = 23 ±4 K. 
A similar value {Trot = 20 K) is obtained with OPR = 3 (fit shown in Fig. [2D and OPR = 2.5 
{Trot = 22 K). The agreement between our model and the Spitzer spectrum is rather good. 
The residuals with respect to the observed spectrum do not show any excess emission in the 
wavelength range 5.8-7.1 fim of water i>2 band emission. 



3.3. Water Production Rates 

From the intensity of the 1^2 H2O band measured on the spectra, we computed the water 
production rates of 71P/Clark and C/2004 Bl (LINEAR) using a Haser model for the water 
density (Tabled]). Because the observations of comet 71P/Clark were affected by pointing 
offsets due to uncertainties in the ephemeris (§ [2]), the production rates retrieved from the 
intensities measured in a small (3.7" x 5.6") and large (3.7" x 17.9") extraction aperture were 
used as a complementary check to the pointing offsets. Consistent values for Q{B.20), equal 
to (1.08 ±0.08) xlO^^ molecules s"^ and (1.04 ±0.08) xlO^^ molecules s"^ for the small and 
large extraction apertures respectively, are obtained applying offsets corresponding to the 
actual position of the nucleus with respect to the targeted position. This derived production 
rate for the 2006 peri helion is consist ent with indirect estimations from the visual brightness 



at previous perihelia (ISekaninalll993l ) 



The water production rate derived for C/2004 Bl (LINEAR) at Vh = 2.21 AU pre- 
perihelion is Q{B.20) = (1.03 ±0.02) xlO^^ molecules s"^ At rn = 2.06 AU post-perihelion, 
the measured Q{}l20) is (4.74 ± 0.42) xlO^'' molecules s~^, which suggests that the comet 
was more active pre-perihelion. This later conclusion is consistent with measurements of the 
total visual magnitude ruy of the comet. For 2005 October 21 and 2006 May 17, i.e., at dates 
close to the Spitzer observations, values m^, = 12.6 and 12.4 are reported, which co rrespond 



to heliocentric visual magnitudes of 11.1 and 11.8, respectively (lGreenll2005l . l2006l ). 



Water production rates were computed assuming a Haser parent-molecule distribution 
for the water density with the parameters given in § 13. 1[ The signal-to-noise ratio achieved 
for C/2004 Bl (LINEAR) pre-perihelion permits the investigation of the radial distribution 
of water throughout the IRS SL2 spectral map (Fig{T^). The intensity of the i>2 band along 
nine extractions is plotted as a function of offset p in Fig. [6](a). The z/2 intensity varies as 
p-0.54 j^Qj-^j^ Qf ^\^Q nucleus which suggests that the water distribution is extended in this 
direction. South of the nucleus, the intensity varies as p~^'^^, which is consistent with a 
Haser distribution. For comparison (Fig. [6^), the surface brightness of the dust emission 
(derived from continuum measurements at 5.45 ± 0.08 /im in each of the nine apertures) 
decreases with increasing offset as p~^-^^ (southern and northern coma averaged). 
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The extended water distribution is best evidenced in the "Q-curve" (cf.. iMumma et al. 



20031 ) shown in Fig. Mjo) which plots the water production required to explain the H2O 



band intensity at the different offsets, assuming a Haser distribution. The apparent Q{B.20) 
increases by more than a factor of 2 from 1 x lO'' km to 4 x 10^ km. This extended 
distribution is suggestive of the presence of subliming icy grains in the coma of C/2004 
Bl (LINEAR). At 2.1 AU from the Sun the lifetime o f dirty ice grains ranges from 10^ — 
10^ s for grains of radii 10^^ — 1 cm (IBeer et al.ll2006l ). This lifetime range corresponds to 
displacements within the coma of typically 10^ — 10^ k m prior to comple te grain sublimation. 
The lifetime of pure ice grains is considerably larger (IBeer et al.ll2006l ). Our measurements 
are most easily explained if pure icy grains or ices mixed with poorly absorbing materials are 
present in the coma of comet C/2004 Bl (LINEAR). However, modeling of such sublimation 
processes is beyond the scope of this paper. 



3.4. Rotational Temperatures 

The measured rotational temperatures (Table [1]) may not correspond to physical gas 
temperatures. As the water molecules expand in the coma, the rotational populations within 
the ground vibrational state relax to a cold fluorescence equilibrium. Thermal equilibrium 
breaks down at distances which depend on the density of the collisional partners (mainly H2O 
molecules and electrons), and which are directly related to the water production rate. Given 
the low water production rate of 71P/Clark and C/2004 Bl (LINEAR), and the observational 
circumstances (field of view, slit offset), we may expect molecules with non-LTE rotational 
populations within Spitzer I RS SL2 slit. Using a mo del of water rotational excitation and 



the methodology explained in lWoodward et al.l (120071 ) . we have computed synthetic 1^2 band 



profiles for the observing conditions of 71P/Clark and C/2004 Bl (LINEAR) from Spitzer 
and the field-of-view corresponding to the extraction apertures. Effective Trot were derived 
by fitting the synthetic spectra with the same procedure as used for the observed spectra. 

For 71P/Clark, the predicted effective Trot at 8" from the nucleus is 11-12 K for a kinetic 
temperature Tkm in the range 10-50 K and Xne = 1, where x„e is a multiplying factor to 
the electron density normalized to the Giotto IP/Halley measurements. If our observations 
had been executed with the sht centered on the nucleus, we would have expected Trot = 25 
to 29 K with these model parameters. These results are for an aperture of 3.7" x 17.9", but 
similar numbers are obtained for a 3.7" x 5.6" aperture: 31-39 K on nucleus, 12-13 K for 
an offset of 8". The size of the collisional region where water molecules are thermalized to 
the gas kinetic temperature is small with respect to the field of view, hence the predicted 
rotational temperature is weakly sensitive to the assumed Tkin- A relatively high Trot is 
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expected for a slit centered on the nucleus because molecules excited by collisions with hot 
electrons have a significant contribution to the signal. Rotational temperatures measured 
in 71P/Clark (14 ± 4 K and 5 ± 13 K for the small and large apertures, respectively) are 
consistent with modeling. 

Model calculations pertaining to the observational circumstances of the C/2004 Bl 
(LINEAR) pre-perihelion data predict Trot = 17-19 K for a 3.7" x 17.9" aperture and kinetic 
temperatures of 10 to 50 K. Note that low gas temperatures may be expected for this weakly 
active (§ 13. 3p comet observed at 2.1-2.2 AU from the Sun. Overall, the agreement with the 
measured Trot value (13.6 ± 2.1 K) is satisfactory. 

In the case of C/2004 Bl (LINEAR) post-perihelion, the expected Trot for a 3.7" x 
17.9" extraction centered on the nucleus is Trot = 17 K assuming T^in = 30 K, with little 
dependency on the assumed T^j^ since water molecules with non-LTE rotational populations 
are sampled. The measured Trot of 23 ± 4 K is consistent with the model within 2-a. A 
better agreement could be obtained by increasing the electron density through the Xne = 1 
parameter. 



3.5. Ortho-to-para Ratio in C/2004 Bl (LINEAR) 



The ortho-to-para ratio of 2.31 ± 0.18 measured in C/2004 Bl (LINEAR) corresponds 
to a spin temperature Tspin = 26^9 K. The ort ho-to-para ratio in water has been mea- 



sured in a dozen of comets (see lBonev et al.ll2007l . for a compilation of measurements before 



2007), and ranges from ~ 1.8 to the equilibrated value of 3 (Tspin from ~ 23 to > 40 K). 
The OPR measured in C/2004 Bl (LINEAR) is in the low ra nge of measured valu es in 
Port cloud comet s , yet comparable to OPR values reported by iMumma et al.l (Il988l ) and 
Dello Russo et al.l J2005h for comets IP/Halley and C/2001 A2 (LINEAR). Our new mea- 
surement for C/2004 Bl (LINEAR) suggests that OPR variation clearly exists among the 
population of Oort cloud comets. 



The meaning of n uclear spin temperatures in cometary molecules is puzzling ( iKawakita et al 



2004J : ICrovisierl 120071 ) . Nuclear spin temperatures may be related to the thermal condi- 
tions experienced by the molecules since their incorporation into cometary nuclei, or con- 
nected to their format ion conditions and subsequent OPR reequilibration in the solar nebula. 
Kawakita et al.l (120041 ) have shown that there is no clear trend between the nuclear spin tem- 
peratures, production rate, heliocentric distanc e and orbital period o f the comet, suggesting 
that the nuclear spin temperature is pristine (IKawakita et al.ll2004l ). The measurement of 
the OPR (and corresponding nuclear spin temperature) in C/2004 Bl (LINEAR) is derived 
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from observations conducted when the comet was at a hehocentric distance, Vh = 2.21 AU. 
This hehocentric distance is intermediate between the range Vh = 0.8-1.2 AU, where most 
OPR measurements were acquired, and the measurement in C/1995 01 (Hale-Bopp) at 2.9 
AU from the Sun. Our Spitzer measurement provides a new data point for future investi- 
gations of Tspin variations with hehocentric distance based on a statistically significant data 
set. 



3.6. Residual Emissions 



The modeled water spectrum wholly accounts for the emission in excess of the dust 
continuum background in the two comets (Figs. [Hand [5]). There is no emission excess at 
6.2 nm that could be attributed to PAH emission. There is not significant evidence for 
moderately broad (AA ~ 0.5 /im) emission peaking at 7.0 /zm (arisi ng from asymmetric 



stretching vibration of C-0) that c ould be due to carbonates minerals (IKemper et al.l 12002 



Crovisier fc Bockelee-Morvan 20081 ). Such emission, if present, would have made the 6.85 fxm 
peak of the H2O band of less contrast than observed. In addition if carbonates where present, 
other narrow (AA ^ 0.2 /xm) emission features near 12.7 fim and 14.8 fim should be present 
jKemper et al.ll2002h . The IRS spectra of comets C/2004 Bl (LINEAR) and 71P/Clark 
over this wavelength interval are featureless. Fig. [71 Thus, carbonates are not present (or 
significantly abundant) in the comet comae created through quiescent sublimation processes 
or dermal jets from isolated active areas on the nucleus based on analysis of Spitzer spectra 
of three comets. 

Excess emission is present at 7.1-7.3 /xm in the spectrum of 71P/Clark, and between 7.15 
and 7.45 fim, peaking at 7.32 /xm in the post-perihelion spectrum of C/2004 Bl (LINEAR) 
(Figs. Hland [5]). Excess emission near 7.3 iim w as also observed in the Spitzer spectrum of 
C/2003 K4 (LINEAR) JWoodward et al.ll2007h . The emission feature observed in C/2004 
Bl (LINEAR) post-perihelion closely resembles that observed for C/2003 K4 (LINEAR), 
but is much more intense relative to the intensity of the water band than in C/2003 K4 
(LINEAR): ~ 10% the intensity of the water band in the post-perihelion spectrum of C/2004 
Bl (LINEAR), ~ 2% in C/2003 K4 (LINEAR) JWoodward et al.ll2007h . Its contrast with 
respect to the continuum background is also more important in C/2004 Bl (LINEAR). 
The feature is not present in the pre-perihelion spectrum of C/2004 Bl (LINEAR). This 
su ggests that the feature is not related to the absolute brightness of the comet. As discussed 



m 



Woodward et al.l (120071 ). the origin of this feature is unclear. Based on examination of 



spectral extractions of standard stars used to cahbrate the IRS campaign, this feature is not 
likely an instrumental artifact, although we cannot completely exclude this possibility. 
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CONCLUSION 



We have presented Spitzer IRS spectra of the water z/2 fundamental vibrational band 
near 6.5 /im in the Jupiter-family comet 71P/Clark and the dynamically new comet C/2004 Bl 
(LINEAR). Modeling of the IR spectra using a fluorescence model of the water emission 
yielded rotational temperature and production rate measurements. The measured low ro- 
tational temperatures are consistent with water excitation models which predict that water 
molecules within the Spitzer IRS slit have partially relaxed (non-LTE) rotational popula- 
tions. The water ortho-to-para ratio (OPR) was constrained to 2.31 ± 0.18 in C/2004 Bl 
(LINEAR) from its pre-perihelion spectrum observed with a high signal-to-noise ratio. The 
corresponding spin temperature Tspin = '26~^l K is in the low range of measured values in other 
Oort cloud comets. In addition, the spatial distributon of water observed (pre-perihelion) in 
the northern coma of C/2004 Bl (LINEAR) is more extended than that expected for nucleus 
outgassing which suggests the presence of subliming icy grains. 

Although the 5.8-7.6 fim spectral domain is rich in spectral signatures of gaseous com- 
pounds and minerals our spectra show no incontrovertible evidence for carbonates, or PAH 
emiss ion. Thus, repor t ed cla ims that these astrobiolog ically significant species are present 



Lisse et al 



2007, 2006 



Bockelee-Morvan et al. 



2007 



remain controversial (i.e., Crovisier fc Bockelee-Morvan 2008 



Woodward et all 120071 ). As seen in comet C/2003 K4 (LIN- 



EAR) (IWoodward et al.l 120071 ). residual emission at 7.1-7.3 fim is present in spectra of 
71P/Clark and C/2004 Bl (LINEAR) whose origin is unclear. 
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1256406, and 1215746 issued by JPL/Caltech to the University of Minnesota. C.E.W. also 
acknowledges support from the National Science Foundation grant AST-0706980. We thank 
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Table 1. WATER IN COMETS: MODEL FITS AND PRODUCTION RATES 



Obs Date 




rh 


A b 

^Spitzer 


Extraction 


Band Intensity'^ 


'^rot 


OPR 






(UT) 


(deg) 


(AU) 


(AU) 


(" X ") 


(10-21 W cni-2) 


(K) 






(10^^ molecules s~^) 


71P/Clark 




















2006 May 27.56 


37.8 


1.57 


0.92 


3.7 X 17.9° 


17.0 ± 1.29 


5 ± 13 


3.0^ 


0.90 


10.43 ± 0.79 


C/2004 Bl (LINEAR) 




















2005 Oct 15.22 


27.5 


2.21 


2.03 


3.7 X 17.9 


13.74 ± 0.28 


13.6 ± 2.1 


2.31 ± 0.18 


0.97 


10.32 ± 0.21 


2006 May 16.22 


28.7 


2.06 


1.61 


3.7 X 17.9 


9.20 ± 0.82 


23 ± 4 


2.3^ 


0.90 


4.74 ± 0.42 



'^Phase angle (defined as vertex angle of Sun-Comet- 6piteer). 

'^Distance to Spitzer spacecraft. 

'^Band intensity between 5.8 /im and 7 /im. 

'^Reduced between 5.85 /im and 7 /im. 

° Offset with respect to the nucleus position by ^ 8" (see §[5]). 

^Assumed. 
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Fig. 1. — Gray-scale surface-brightness images derived from Spitzer IRS spectral mapping 
data cubes of comets C/2004 Bl (LINEAR) and 71P/Clark: (a) IRS short-low order 2 (SL2 
5.3-7.5 /im) observation of comet C/2004 Bl (LINEAR) at Vh = 2.2 AU (pre-perihelion) 

(b) IRS SL2 observation of comet C/2004 Bl (LINEAR) at r?, = 2.1 AU (post-perihelion) 

(c) IRS SL2 observation of comet 71P/Clark at = 1.6 AU (pre-perihelion); (d) IRS long- 
low order 2 (LL2; 14-20 /im) observation of comet 71P/Clark obtained as part of the same 
Spitzer AOR used to acquire the SL2 observation of panel (c). Arrows in panels (a), (b), and 
(c) delineate the projected sun (©) and comet velocity (v) directions. The 71P/Clark SL2 
observation did not cover the comet nucleus; however, the LL2 observation did. To emphasize 
this point, we have marked the JPL/Horizons ephemeris position of comet 71P/Clark with 
a circle, and placed outlines of the LL2 and SL2 modules in panels (c) and (d), respectively. 
Spectral extractions discussed in § |2] are centered on the peak of the surface brightness 
distribution. 



- 19 - 




Fig. 2.— Model fits to the spectra of 71P/Clark (A, B) and C/2004 Bl (LINEAR) pre- 
perilielion (C, D) and post-perilielion (E, F) observed witli Spitzer on 2006 May 27.56 UT, 
2005 October 27.56 UT and 2006 May 16.22 UT, respectively. The extraction aperture 
is 3.7" X 17.9". Data are shown in black, with the model fits of water emission (in red) 
superimposed. For 71P/Clark and C/2004 Bl (LINEAR) post-perihelion, model fitting was 
performed with the rotational temperature T^ot taken as a free parameter, and the ortho- 
to-para ratio set to OPR = 3. In contrast, for the pre-perihelion spectrum of C/2004 Bl 
(LINEAR), the fit was obtained with both the OPR and Trot taken as a free parameters. The 
underlying continuum, described by a polynomial of degree 5, was also fitted simultaneously. 
In panels B and D, the fitted continuum background has been subtracted. 
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Fig. 3. — Contours showing the variation of with Trot and OPR in the region that 
provides the best fit to the C/2004 Bl (LINEAR) pre-perihehon spectrum. Iso-x^ contour 
levels correspond to increases in by steps of 1 from the minimum value (denoted by small 
black dot), (a) The contours are calculated by holding the polynomial coefficients and the 
water production rate fixed at their optimum values, (b) The contours shown in red color 
are calculated by allowing those parameters to vary to minimize fo^' each pair of values 
of Trot and OPR. Black contours are those shown in plot (a). 
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Fig. 4. — Model fit to the 2006 May 27.56 UT (pre-perihelion) spectra of comet 71P/Clark 
with extraction apertures of 3.7"x 5.6" (top) and 3.7"x 17.9" (bottom). Data are shown in 
black dots with errorbars. The synthetic water spectrum obtained from model fitting with 
the ortho-to-para ratio set to OPR = 3 is shown in red. The retrieved Trot is 14 ± 4 K and 
5 ± 13 K for the small and large apertures, respectively. The residual spectrum is shown 
in red on the bottom, on which are su perimposed iii arbitrary units an interstellar PAH 
spectrum typical of class A sources frorn iPeeters et al.l (120021 ) {blue spectrum), and a model 
of carbonate emission from Lisse et all J2006h {green spectrum). 
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Fig. 5.— Model fits to the 2005 October 27.56 UT (pre-perihelion, top) and 2006 May 
16.22 UT (post-perihelion, bottom) spectra of comet C/2004 Bl (LINEAR) obtained with 
an extraction aperture of 3.7" x 17.9". Data are shown in black dots with errorbars, while 
other colored curves are as described in Fig. |H Model fitting to the pre-perihelion spectrum 
yields Trot = 13.6 ± 2.1 K and OPR = 2.31 ± 0.18. The best fit obtained with the OPR 
fixed to OPR = 3 {Trot = 13.4 K) is shown by a black dotted line. For the post-perihelion 
spectrum, the OPR is fixed to a value of 2.3, and the derived rotational temperature is 
Trot = 23 ± 4 K. 
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Fig. 6. — Distribution of water and dust as a function of offset, p, with respect to the 
nucleus position in C/2004 Bl (LINEAR) at Vh = 2.2 AU (pre-perihehon) . Negative offsets 
correspond to southern extractions with respect to the peak in the surface brightness of 
the coma, Fig. [T](a). (a) Total water z/2 band intensity between 5.8 and 7 /im {black circles) 
and dust continuum {red circles) surface brightness determined from equivalent narrow-band 
continuum (divided by a factor of 17) measured at A = 5.45 ±0.08 /im. Extraction apertures 
are 5.6" x 5.6" , except for the most central extraction which has an 5.6" x 11.1" aperture. 
For clarity, the central (offset = —1") dust point has been shifted to an offset equal to 
zero. The water intensity varies as p^^'^^ for the southern extractions, and as p~^-^^ for the 
northern extractions, while the dust emission (average of the northern and southern coma 
values) decreases with increasing offset as p~^-^°. (b) The Q-curve, i.e., the apparent Q{B.20) 
production rates required to explain the H2O band intensity at the different offsets in the 
assumption of a Haser parent molecule distribution. 
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Fig. 7.— The Spitzer IRS short-high (SH) spectra of comet C/2004 Bl (LINEAR) over 
the wavelength range from 12.0 fim to 16.0 fim where narrow resonances from the out- 
of-plane and in-plari e bending mode of the (COs)^^ ion in carbonate species are expected 
( iKemper et al.ll2002l ). The grey vertically shaded areas indicate the wavelength locations of 
the spectral order overlap in the instrument and the discontinuity of the spectra in these 
regions are artifacts. The IRS spectra are featureless in the wavelength regions of the car- 
bonate resonances (narrow features near 12.7 fim and 14.8 /im), suggesting that carbonates 
are not present with large abundance, if at all, in the emitting coma materials. 



